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ABSTRACT 

We discuss a new method for inferring the stellar mass of a distant galaxy of known redshift 
based on the combination of a near-infrared luminosity and multi-band optical photometry. The 
typical uncertainty for field galaxies with I < 22 in the redshift range < z < 1 is a factor 
2. We apply this method to a newly-constructed sample of 321 field galaxies with redshifts and 
Hubble Space Telescope morphologies enabling us to construct the stellar mass density associated 
with various morphologies as a function of redshift. We find a marked decline with time in the 
stellar mass associated with peculiar galaxies accompanied by a modest rise in that observed for 
ellipticals. The result suggests that peculiars decline in abundance because they transform and 
merge into regular systems. The star formation rate per unit stellar mass indicates that massive 
systems completed the bulk of their star formation before redshift one, whereas dwarfs continue 
to undergo major episodes of activity until the present epoch. 

Subject headings: galaxies: evolution, galaxies: fundamental parameters, galaxies: stellar content 



1. Introduction 

Our understanding of the evolution of galaxies 
during the important redshift range < z < 1 
has advanced significantly in recent years through 
the completion of statistically-complete redshift 
surveys (Lilly et al 1995, Ellis et al 1996) and 
the association of these data with morphological 
information from Hubble Space Telescope (HST, 
Brinchmann et al 1998, B98). This progress has 
been mirrored by the formulation of detailed semi- 
analytic models which attempt to reproduce the 
observational trends in the context of hierarchical 
galaxy formation (Baugh et al 1998; Kauffmann 
et al 1999). 

One of the most intriguing results from HST 
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imaging and ground-based spectroscopy is that the 
steep slope in the blue number counts arises from 
a surprisingly abundant population of irregular, 
star-forming galaxies whose local counterparts are 
unclear (Glazebrook et al. 1995). Such galaxies 
contribute a negligible amount to the local rest- 
frame blue luminosity density (Marzke et al 1998) 
but equal the contribution of regular spirals at red- 
shifts z ~1 (B98). In hierarchical models, "mor- 
phologies" transform according to the mass char- 
acteristics of mergers. Recent merging leads to 
growth in regular systems at the expense of lower 
mass irregulars (Baugh et al 1998). 

It is observationally difficult to verify this 
hypothesis using traditional methods based on 
magnitude-limited morphological counts, not only 
because galaxies transform from one kind to an- 
other but also since they presumably brighten 
in optical luminosity as they merge. To make 
progress we require an independent 'accounting 



variable' capable of tracking the likely assembly 
and transformation of galaxies during the interval 
< z < 1. The colour and emission line charac- 
teristics are transient properties and poorly-suited 
for this purpose (Glazebrook et al 1999). Ideally 
we seek a quantity that is readily observable but 
unaffected by morphological transformations and 
mergers. 

The dynamical or stellar mass is the obvious 
choice. The former is directly connected to theo- 
retical predictions but difficult to observe except in 
favourably oriented regular systems (e.g. Vogt et 
al 1997). The stellar mass is less well-determined 
by theory but has the very considerable advan- 
tage that it can be derived robustly for complete 
samples of galaxies arranged morphologically us- 
ing multi-band optical and near-infrared (near-IR) 
photometry enabling a complete inventory to be 
made for each type at various epochs. Since stel- 
lar mass is either conserved, or a slowly increas- 
ing quantity, it enables us to track the fate of the 
abundant z ~ 1 irregular systems. 

This paper is concerned with demonstrating the 
promise of this technique and its initial application 
to a statistically-complete sample of field galaxies 
of known redshift. We have adopted a cosmologi- 
cal model with f2 = 0.35, Ao = 0.65 and h = 0.65 
throughout but our results are not sensitive to this 
choice. 

2. Infrared Photometry for a Complete 
HST-based Sample 

Our dataset is an extension of that introduced 
by B98 which resulted from HST imaging of the 
CFRS and LDSS redshift surveys. The key fea- 
tures are spectroscopic redshifts for a statistically- 
complete sample with WFPC-2 7 8 i4 imaging suit- 
able for morphological classifications, both visu- 
ally and with automated classifiers (Abraham et 
al 1996). By adding archival data for the Hubble 
Deep Field North (HDF-N) and its flanking fields 
(Williams et al 1996) and the Hawaii Deep Fields 
(Cowie et al 1996), we obtain a spectroscopic and 
morphological sample of 498 field galaxies with 
17.5 < Iab < 22.5. Further details of this sample 
and its uniformity are discussed by Brinchmann 
(1999). 

Here we analyse a subset of 321 galaxies for 
which deep K-band (2.2/xm) photometry has been 



obtained from the 3.8m UK Infrared Telescope 
and the literature. Imaging was performed by the 
authors using IRCAM3 during three photometric 
nights 21-24/5 1998 using a K filter. Additional 
data came from Menanteau et al's (1999) HK' sur- 
vey of the Medium Deep Fields, Dickinson's (un- 
published) K exposures of HDF-N and Barger et 
al's (1999) HK' imaging of the HDF flanking fields 
and some original CFRS K' photometry (Hammer 
et al. 1997). Total magnitudes & aperture colours 
have been re-measured in a consistent manne; the 
various IR filter differences are taken into account 
in the SED fitting algorithm discussed below. 

3. An Improved Estimator of Stellar Mass 
for Distant Galaxies 

Many authors have pointed out that the in- 
frared luminosity of a galaxy is a more robust es- 
timator of its stellar mass than that derived at 
optical wavelengths. The key role of if-band mag- 
nitudes was emphasised as a possible test of the 
merging hypothesis by Broadhurst, Ellis & Glaze- 
brook (1992) and this has motivated numerous 
infrared-based deep surveys (Cowie et al 1996, Co- 
hen et al 1999). The K luminosity distribution of 
distant galaxies was discussed in support of the hi- 
erarchical picture by Kauffmann & Chariot (1998) 
who demonstrated its insensitivity to the previous 
star formation history. The main uncertainty in 
the inferred stellar masses arises from the age of 
the stellar population (Rix & Rieke 1993). 

The basis of our new method is to apply a 
correction for this uncertainty. We adopt a grid 
of evolutionary synthesis models with simple star 
formation histories (Bruzual and Chariot 2000, 
BC, Table 1) and select the most appropriate evo- 
lutionary history by optimal fitting of the opti- 
cal+IR photometry. Extinction, age and metallic- 
ity are allowed to vary within the bounds tabu- 
lated in a manner similar to that adopted for pho- 
tometric redshift determination (c.f. Giallongo et 
al 1998) and the mass is estimated by normalising 
the best fit SED to the observed -ftf-band magni- 
tude. 

Clearly there are some significant assumptions 
in this approach. Our approach here is to intro- 
duce the method and its potential whilst stressing 
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Table 1 

Grid of Model Spectral Energy Distributions 



Parameter Range 

Age 10 8 -2 x 10 10 years in steps of 0.1 in log age. 

SFR r =0.1, 0.5, 1.0, 2.5, 5.0 Gyr and constant 

Metallicity Z = 0.02, Z = 0.004 

E(B - V) 0.0, 0.2, 0.4 — Calzetti (1997) law 



12 r 




Log M(dyn) 



Fig. 1. — A comparison of dynamical mass estimates within three scale-radii (Md yn = V 2 R/G) with the 
stellar mass estimates detailed in the text. Symbols with error bars: six galaxies from Vogt et al. (1997) from 
the HDF flanking fields; these have magnitudes measured on a uniform system. Filled circles: 30 galaxies 
from Boselli et al. (1997). The inset histogram shows the average mass distribution from the Monte Carlo 
simulations described in the text and compares this with a Gaussian of width 0.3 in logM star . 



3 



that independent verification is clearly desirable. 
We assume (a) that the IMF does not vary sys- 
tematically with morphological type, and b) that 
the BC code gives reasonable eastimates of the 
near-IR fluxes for a range of stellar populations. 
Both assumptions are, to some extent, supported 
by data for external galaxies (e.g. Kennicutt 1998) 
and comparisons of spectral synthesis codes by 
Chariot, Worthey & Bressan (1996). For much of 
the discussion, relative stellar masses will suffice 
but where absolute masses and estimates of star 
formation rates are required we adopt a Salpeter 
IMF with lower and upper mass cut-offs of 0.1 M 
and 125 M . 

To assess the random errors arising from our 
method, we propagate magnitude errors through 
the fitting procedure using a Monte Carlo method. 
This gives an intrinsic error on the mass estimates 
with a median error of 0.20dex in log M st ar (see in- 
set histogram in Figure 1). Further errors might 
arise due to our choice of exponentially declining 
star formation histories, and to quantify this we 
have explored the sensitivity of the method to a 
secondary burst of SF and find that this can in- 
troduce scatter of ~ 0.2 in logM sta r with no sign 
of systematic offsets. Taken together this gives an 
estimate of the theoretical uncertainty in the mass 
estimates of « 0.3dex in logM sta r- 

One potentially useful test of our method 
is a comparison of our stellar mass estimates 
with those determined dynamically using rotation 
curves. We currently have a small overlap of 6 
galaxies with Vogt et al (1997) and assuming that 
these have accurately estimated terminal velocity 
speeds we can estimate the dynamical mass within 
three scale-radii using Md yn = V 2 R/G. The com- 
parison of this mass estimate with the stellar mass 
is shown in Figure 1. The same exercise has been 
performed with the local sample of Boselli et al 
(1997) although we expect a larger scatter here as 
the photometric scales have not been accurately 
aligned. Figure 1 shows both masses correlate 
with a scatter consistent with earlier estimes. Of 
course, the relationship between dynamical and 
stellar mass depends on the physical properties 
and is of intrinsic interest. The mean offset cor- 
responds to M star — 0.3-0.5 Mdyn- Clearly it will 
be instructive to construct Figure 1 for a larger 
sample. 



4. Mass Assembly and Specific Star For- 
mation Rates 

Although considerable progress has been made 
in delineating the star formation history of field 
galaxies (e.g. Madau, Pozzetti & Dickinson 1998), 
such probes remain uncertain particularly as there 
is not yet good agreement on what represents a 
reliable diagnostic. Moreover, it is clear from 
the considerable diversity of predictions in semi- 
analytic models (c.f. Baugh et al 1998) that cos- 
mic star formation histories cannot precisely test 
hierarchical theories of galaxy formation. Whilst 
considerable effort has been invested in incorpo- 
rating realistic aspects of star formation and feed- 
back into numerical models of galaxy evolution, 
here we explore an alternative route, viz. using 
our method to bring the empirical data closer to 
the fundamental aspects of hierarchical theory by 
constructing the redshift dependence of the stellar 
mass function. 

To complement the co-moving volume-averaged 
star formation history p*(z), we seek to construct 
its stellar mass equivalent p m (z) and particularly 
its dependence on galaxy morphology. With suffi- 
cient data we could also construct the stellar mass 
function $(m), analogous to the luminosity func- 
tion, using a 1/V max formalism. Its consideration 
allows us to discuss the important question of pos- 
sible incompleteness which would otherwise under- 
estimate the stellar mass density p m (z). 

Incompleteness can arise in various ways. The 
most critical is that galaxies beyond our I&u limit 
in a given redshift range contribute to the infrared 
flux and hence the stellar mass density. Ultimately 
we seek a if -limited version of our current sam- 
ple from which galaxies can, to first order, be se- 
lected to a fixed stellar mass limit at a given red- 
shift. Meanwhile, to assess this incompleteness, 
we utilise the form of the K-band galaxy lumi- 
nosity function (LF) obtained by combining the 
recent infrared survey of Cohen et al (1999) with 
the if-band LF of Loveday (1999). 

Assuming a fixed mass/light ratio of M stal /LK=0.8, 
this LF suggests that if we restrict ourselves to 
galaxies with stellar masses in the range 10.5 < 
logM s tar < 11.6, we will be > 80% complete in p m 
across the interval 0.2 < z < 1.0, with the great- 
est incompleteness in the highest redshift bin. In 
the subsequent discussion we have restricted our 
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Fig. 2. — The integrated stellar mass density for galaxies with 10.5 < logM sta r < 11-6 as a function of 
redshift and visual morphology. The shaded regions show the predictions of the mass density in peculiar and 
elliptical galaxies from the simple merging models described in the text. 



analysis to this mass range and taken the incom- 
pleteness corrections and their uncertainty into 
account. 

Figure 2 shows the comoving stellar mass den- 
sity in galaxies partitioned by morphology as a 
function of redshift with error bars derived by 
bootstrap resampling methods. The solid horizon- 
tal line shows the total local stellar mass density 
inferred from a recent review by Fukugita, Hogan 
& Peebles (1998) adjusted to our IMF. The most 
significant feature is the remarkable decline in 
mass density associated with galaxies of irregular 
morphology. This is complemented by a modest 
increase in that associated with regular galaxies, 
particularly spheroidals, ensuring that the mean 
stellar mass density remains approximately con- 
stant at all redshifts probed. Inspection of the 
mass function for high redshift irregulars (Brinch- 
mann 1999) shows a significant number harbour 
stellar masses as large as those for the most mas- 
sive regular galaxies suggesting that these simply 
transform into the latter systems perhaps as a re- 
sult of the onset of more regular star formation 
patterns. However, the decline in the irregular 



galaxy luminosity density may also be driven by 
merging, causing a gradual increase in the mass 
density associated with regular galaxies. 

What rate of merging is necessary to explain 
Figure 2? By accounting for evolution in terms 
of stellar mass, rather than blue luminosity, we 
can consider this issue using recent estimates of 
the redshift-dependent merger fraction (Le Fevre 
et al 2000, and references therein). Adopting a 
merger fraction oc (1 + z) 3A , Brinchmann (1999) 
finds a good fit to the morphologically-dependent 
redshift distribution and luminosity functions if 5- 
10% of the present population of E/SOs formed 
via mergers over < z < 1. This solution and 
its uncertainty is illustrated by the shaded area in 
Figure 2 where the upper locus refers to a model 
where the merger progenitors are assumed to be 
of average mass and with colours similar to local 
Sbc galaxies over the entire redshift range. How- 
ever the data indicate that peculiar galaxies be- 
come less massive at low redshift, hence we show 
a model where merger progenitors at z = 0.5 are 
30% less massive than those at z = 1 as the lower 
locus. The total mass density is normalised to 



5 



-9 
-10 
-11 



fa 
m 



3 -12 



-13 



-14 



— "Doubling " 



Median error 



4 . ^TV^/ 



—"10% Increase" 



• 0.2<Z<0.50: 
. 0.5<Z<0.75: 
■ 0.75<Z<1.0: 



1.70 -1.19 Log M 
1.96 -1.16 Log M 

1.71 -1.12 Log M 



0.1 
1 

10 



O 

100^ 



1000 



10" 



i 1CT 



8 9 10 11 12 13 

T*og M lurn 

Fig. 3. — The specific star formation rate, R, for the galaxies in the sample. The three lines show orthogonal 
least-squares fit to the three redshift ranges indicated. The fit parameters are shown in the lower left corner. 
The right hand panel shows the doubling time in Gyr, assuming constant star formation. Arrows show 2a 
upper-limits for galaxies with no detected [On] emission. 



the observed amounts in each redshift range, and 
the model chosen such that it reproduces the lo- 
cal morphologically-split LFs from Marzke et al. 
(1998). The increase with time in the E/SO mass 
density is well matched in both models whereas 
the decrease in the peculiar mass fraction is over- 
estimated. Converting some peculiar galaxies into 
spirals (whose mass density was assumed to be 
constant in the foregoing) would improve the fit. 

A further important result can be extracted 
from Figure 2, namely that the total stellar mass 
density in massive galaxies is approximately con- 
stant over the interval 0.2 < z < 1. This, in 
turn, implies that most of the stellar content of 
these galaxies must have formed at higher redshift. 
At first sight this seems a surprising result given, 
over the same redshift range, the volume-averaged 
star formation rate remains high (Madau et al. 
1998). Our results can be more easily understood 
by considering the contribution that star forma- 
tion makes to the total mass of a galaxy. 

Consider the specific star formation rate R = 
SFR/M stal which defines a useful measure of the 
rate at which new stars add to the assembled mass 



of a galaxy. For massive galaxies, a much larger R 
is needed to make a fractionally significant mass 
contribution than would be the case for a low mass 
system. Figure 3 shows R as a function of stellar 
mass for the CFRS galaxies in our survey (Ham- 
mer et al. 1997). Different symbols relate to three 
redshift ranges, arrows show 2a upper limits. Star 
formation rates were derived using the luminosity 
observed in the [O II]3727 A emission line using 
(Guzman et al. 1997) 

SFRiMoyr- 1 ) = i -^a-oa(m b -m b , q ) EW ([On])^ 

^ (1) 

where Mb is the blue luminosity. Note that Fig- 
ure 3 makes a strong statement about the upper 
envelope of the distribution where incompleteness 
has no effect. 

The derivation of stellar masses for each galaxy 
in our survey enables us to ask when galaxies of a 
particular mass range formed most of their stars. 
The right hand axis in Figure 3 shows the dou- 
bling time for the stellar mass of a galaxy assum- 
ing constant star formation, with the horisontal 
lines indicating the value of R required to double 
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and increase by 10% the stellar mass of a galaxy 
from z = 1 to the present. Clearly massive galax- 
ies for the bulk of their stellar content at z > 1. 
At lower redshift, the pattern continues for lower 
mass galaxies consistent with the "downsizing" 
picture first introduced by Cowie et al (1996) and 
our findings are also in very good agreement with 
Guzman et al. (1997) when the differences in mass 
estimation are taken into account. Given our in- 
completeness, dwarf galaxies may well form over 
a wide redshift range since we cannot probe such 
mass ranges at all redshifts. 

A potential problem with this argument would 
be massive amounts of obscured star formation, 
since we have made no dust correction beyond the 
local calibration. We cannot rule this out conclu- 
sively, but observe that only few galaxies in our 
sample show significant far-IR emission (Flores et 
al. 1999), and even a constant star formation rate 
of 10M Q /yr will only cause a 10% increase in the 
mass density over the redshift and mass range we 
here consider. 

5. Conclusions 

We have presented a new method for deter- 
mining reliable stellar masses for field galaxies of 
known redshift and applied it to a uniform sam- 
ple of 321 galaxies with HST morphologies. We 
demonstrate via comparisons and simulations that 
our method yields stellar masses with a precision 
of logAM s tar ~ 0.3dex, although further calibra- 
tions and comparisons with dynamical estimates 
would be valuable. The main advantage of our 
method is that, with few assumptions, it provides 
stellar masses with greatly reduced observational 
requirements and with no bias against galaxy type 
or orientation. It is therefore particularly well- 
suited for considering the mass density of galaxies 
of different types as a function of redshift. 

Whereas only a small proportion (< 10%) of 
massive field spheroidal galaxies can have formed 
by merging between z = and z = 1, this most 
likely arises from a significant demise in the mass 
density of high-z peculiar galaxies. The transfor- 
mation of irregulars into more massive regular sys- 
tems is shown to be in reasonable agreement with 
the evolution of the merger fraction found inde- 
pendently from pair statistics. 

We find no evidence for any change in the mass 



density of spiral galaxies, suggesting the merg- 
ing of spiral galaxies into elliptical galaxies be- 
tween z = 1 and the present cannot be a dom- 
inant process. Moreover, the total stellar mass 
density hardly increases over the range sampled. 
Both results support the contention that most 
massive galaxies formed their stars prior to or 
around z = 1. Lower mass galaxies remains ac- 
tive to much more recent epochs consistent with 
the " downsizing" picture introduced by Cowie and 
collaborators. 
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